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Abstract 

Valiant introduced matchgate computation and holographic algorithms. A number of seemingly exponential 
time problems can be solved by this novel algorithmic paradigm in polynomial time. We show that, in a 
very strong sense, matchgate computations and holographic algorithms based on them provide a universal 
methodology to a broad class of counting problems studied in statistical physics community for decades. They 
capture precisely those problems which are #P-hard on general graphs but computable in polynomial time on 
planar graphs. 

More precisely, we prove complexity dichotomy theorems in the framework of counting CSP problems. 
The local constraint functions take Boolean inputs, and can be arbitrary real- valued symmetric functions. We 
prove that, every problem in this class belongs to precisely three categories: (1) those which are tractable 
(i.e., polynomial time computable) on general graphs, or (2) those which are #P-hard on general graphs but 
tractable on planar graphs, or (3) those which are #P-hard even on planar graphs. The classification criteria 
are explicit. Moreover, problems in category (2) are tractable on planar graphs precisely by holographic 
algorithms with matchgates. 

1 Introduction 

Given a set of functions J 7 , the Counting Constraint Satisfaction Problem #CSP(J r ) is the following problem: An 
input instance consists of a set of variables X — {x\, X2, ■ ■ ■ , x n } and a set of constraints where each constraint is a 
function / € T applied to some variables in X. The output is the sum, over all assignments to A, of the products 
of these function evaluations. This sum-of-product evaluation is called the partition junction. In the special case 
where / e T outputs values in {0, 1} it counts the number of satisfying assignments. But constraint functions 
taking real or complex values are also interesting, called (real or complex) weighted #CSP. Our T consists of 
real or complex valued functions in general. There is a deeper reason for allowing this generality: The theory of 
holographic reductions is a powerful tool which operates naturally over C, even if the original problem has only 
0-1 valued functions. 

A closely related framework for locally constrained counting problems is called Holant Problems [2J [H] . This 
framework is inspired by the introduction of Holographic Algorithms by L. Valiant 40, 39J. In two ground-breaking 
papers [351 HO] Valiant introduced matchgates and holographic algorithms based on matchgates to solve a number 
of problems in polynomial time, which appear to require exponential time. At the heart of these exotic algorithms 
is a tensor transformation from a given problem to the problem of counting (complex) weighted perfect matchings 
over planar graphs. The latter problem has a remarkable P-time algorithm (FKT- algorithm) [33, 24, 25 . Planarity 
is crucial, as counting perfect matchings over general graphs is #P-hard [35] . Most of these holographic algorithms 
use a suitable linear basis to realize locally a symmetric function with at most 3 Boolean variables on a matchgate. 
This work has been extended in [11 . In particular we have obtained a complete characterization of all realizable 
symmetric functions by matchgates over the complex field C. 

The study of "tractable #CSP" type problems has a much longer history in the statistical physics community 
(under different names). Ever since Wilhelm Lenz who invented what is now known as the Ising model, and 
asked his student Ernst Ising [33] to work on it, physicists have studied so-called "Exactly Solved Models" [TII3U]. 
In the language of modern complexity theory, physicists' notion of an "Exactly Solvable" system corresponds 
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to systems with polynomial time computable partition functions. This is captured completely by the computer 
science notion of "tractable #CSP" . In Physics, many great researchers worked to build this intellectual edifice, 
with remarkable contributions by Ising, Onsager, C.N.Yang, T.D.Lee, Fisher, Temperley, Kasteleyn, Baxter, 
Lieb, Wilson etc [121 Ell H31 EH [2H E31 12H ESI III UHl IH] ■ A central question is to identify what "systems" 
can be solved "exactly" and what "systems" are "difficult". The basic conclusion from physicists is that some 
"systems" , including the Ising model, are "exactly solvable" for planar graphs, but they appear difficult for higher 
dimensions. There does not exist any rigorous or provable classification. This is partly because the notion of 
a "difficult" partition function had no rigorous definition in physics. However, in the language of complexity 
theory, it is natural to consider the classification problem. In this paper we do that, in the more general setting 
of #CSP with real valued constraint functions. This will also shed light on why the valiant efforts by physicists 
to generalize the "exactly solved" planar system to higher dimensions failed. (In the appendix we will give some 
more background.) 

Now turning from Physics to CS proper, after Valiant introduced his holographic algorithms with matchgates, 
the following question can be raised: Do these novel algorithms capture all P-time tractable counting problems on 
planar graphs, or are there other more exotic algorithmic paradigms yet undiscovered? A suspicion (and perhaps 
an audacious proposition) is that they have indeed captured all tractable planar counting problems. If so it would 
provide a universal methodology to a broad class of counting problems studied in statistical physics and beyond. 
The results of this paper can be viewed as an affirmation of that suspicion. Within the framework of weighted 
Boolean ^CSP problems our answer is YES, for all symmetric real valued functions. 

While #CSP problems provide a natural framework to address this question, it turns out that the deeper 
reason comes from Holant problems, which can be described as follows: An input graph G = (V, E) is given, 
where each v g V is attached a function f v e J 7 , mapping {0, l} de s( 1 ') — ► R or C. We consider all edge assignments 
a : E — > {0, 1}. For each a, f v takes its input bits from the incident edges E(v) at v, and evaluates to f v (a \e{v))- 
The counting problem on instance G is to compute Holantg = X)a L~Ley ■f v ( <J l-EO))- f n effect, in a Holant 
problem, edges are variables and vertices represent constraint functions. This framework is very natural; e.g., the 
problem of Perfect Matching corresponds to attaching the Exact-One function at each vertex, taking 0-1 
inputs. The class of all Holant problems with function set F is denoted by Holant(.F). 

Every #CSP problem can be simulated by a Holant problem. Represent any instance of a #CSP problem 
by a bipartite graph where LHS are labeled by variables and RHS are labeled by constraints. Denote by =&: 
{0, l} k — > {0, 1} the Equality function of arity k, which is 1 on fe and l k , and is elsewhere. Then we can 
turn the #CSP instance to an input graph of a Holant problem, by replacing every variable vertex v on LHS by 
=deg(u)- In fact, #CSP(J") is exactly the same as Holant (J 7 U {=*| k > I}). Thus, #CSP problems can be viewed 
as Holant problems where all Equality functions are available for free, or assumed to be present. However, 
when we wish to discuss some restricted classes of counting problems, e.g., for 3-rcgular graphs, the framework of 
Holant problems is the more natural one. And as it turns out, the main technical breakthrough for our dichotomy 
theorem for planar #CSP comes from Holant problems. 

In this paper we will only consider Boolean variables X. For a symmetric function on k variables, we denote 
it as [/o, /i, . . . , fk], where is the value of / on inputs of Hamming weight i. E.g., (=i) = [I, I], (= 2 ) = [I, 0, 1] 
and (=3) = [1, 0, 0, 1] etc. When we relax Holant problems by allowing all Equality functions for free, we obtain 
#CSP. We can also consider other relaxations. Let = [1,0] and 1 = [0, 1] denote the constant and 1 unary 
(arity 1) functions. Then Holant c is the natural class of Holant problems where and 1 are free. This amounts 
to computing Holant on input graphs where we can set or 1 to some dangling edges (one end has degree 1). 
Another class of Holant problems is called Holant* problems where we assume all unary functions [ito,ui] are 
free. 

In [IB] we obtained a dichotomy theorem for (complex) Holant* problems and (real) Holant c problems. The 
dichotomy criterion for Holant* problems is still valid for planar graphs. The proof of dichotomy theorems in this 
paper starts from there. 

In Section |4j we prove that for any real- valued symmetric function set J-, the planar Holant c (J r ) problem 
is tractable (i.e., computable in P) iff cither it is tractable over general graphs (for which we already have an 
effective dichotomy theorem [16] ). or it is tractable because every function in J- is realizable by a matchgate, in 
which case the planar Holant c (J r ) problem is computable by matchgates in P-time using FKT. In all other cases 
the problem is ^P-hardQ A crucial ingredient of the proof is a crossover construction whose validity is proved 

* Strictly speaking, we must only consider T where functions take computable real numbers; this will be assumed implicitly. 
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algebraically, which seems to defy any direct combinatorial justification. 

Our second theorem (Section [5]) is about planar #CSP problems. We prove that for any set of real-valued 
symmetric functions J 7 , the planar #CSP(J r ) problem is tractable iff either it is tractable as #CSP(J r ) without 
the planarity restriction (for which we have an effective dichotomy theorem |16j). or it is tractable because 
every function in T is realizable by a matchgate under a specific holographic basis transformation. Thus planar 
#CSP(J r ) is solvable by a holographic algorithm in the second case. For all other T the problem is #P-hard. 
The proof of this dichotomy theorem for planar #CSP is built on the one for planar Holant c in SectionH] 



Our third result is a dichotomy theorem for planar 2-3 regular bipartite Holant problems (Section |6|T (This 
theorem deals with Holant problems without assuming unary and 1.) This includes Holant problems for 3- 
regular graphs as a special case. The tractability criterion is the same: Either it is tractable for general graphs 
(for which we also have an effective dichotomy theorem [9!), or it is tractable by a suitable holographic algorithm, 
which is a holographic reduction to FKT using matchgates. In all other cases the problem is #P-hard. 

The three dichotomy theorems are not mutually subsumed by each other and are of independent interest. In 
each framework the respective theorem is a demonstration that holographic algorithms with matchgates capture 
precisely those ^P-hard problems which become tractable for planar graphs. 

2 Preliminaries 

2.1 Problem and Definitions 

Our functions take values in C by default. The framework of Holant problems is defined for functions mapping 
any [q] k — > C for a finite q. Our results in this paper are for the Boolean case q = 2. So we give the following 
definitions only for q = 2 for notational simplicity. 

A signature grid Q = (H, J 7 , 7r) consists of a graph H = (V, E), and a labeling ix which labels each vertex with 
a function f v e T . The Holant problem on instance f2 is to compute Holantsi = Tiv^v fv(°~ |e(u))> a sum over 
all edge assignments a : E — > {0, 1}. A function f v can be represented as a vector of length 2 dcg ( v \ or a tensor 
in (C 2 )® dcg ( u '. A function / e T is also called a signature. We denote by the Equality signature of arity k. 
A symmetric function f on k Boolean variables can be expressed by [fo, fi, ■ ■ • , fk], where fa is the value of / on 
inputs of Hamming weight i. Thus, (=k) — [1, 0, . . . , 0, 1] (with k — 1 zeros). A Holant problem is parameterized 
by a set of signatures. 

Definition 2.1. Given a set of signatures J-, we define a counting problem Holant(J-~); 
Input: A signature grid Q — (G,J-,ir); 
Output: Holant^. 

Planar Holant problems are Holant problems on planar graphs. 

Definition 2.2. Given a set of signatures J 7 , we define a counting problem Pl-Holant(J r ): 
Input: A signature grid Q = (G,J-,tt), where G is a planar graph; 
Output: Holanto- 

We would like to characterize the complexity of Holant problems in terms of its signature sets. [^] For some 
J 7 , it is possible that Holant(J r ) is #P-hard, while Pl-Holant(J r ) is tractable. These new tractable cases make 
dichotomies for planar Holant problems more challenging. This is also the focus of this work. Some special 
families of Holant problems have already been widely studied. For example, if T contains all Equality signatures 
{=1, =2, =3, • ■ •}, then this is exactly the weighted #CSP problem. P1-#CSP denotes the restriction of #CSP 
to planar structures, i.e., the standard bipartite graphs representing the input instances of #CSP are planar. In 
[16j , we also introduced the following two special families of Holant problems by assuming some signatures are 
freely available. 

Definition 2.3. Let U denote the set of all unary signatures. Given a set of signatures J- , we use Holant* (J 7 ) 
(or Pl-Holant*( T) respectively) to denote Holant {T U U) (or Pl-Holant( T U U) respectively). 

^Usually our set of signatures J 7 is a finite set, and the assertion of either Holant(J 7 ) is tractable or #P-hard has the usual 
meaning. However our dichotomy theorem is actually stronger: we allow T to be infinite, e.g., to include {=i, =2, =3, ■ ■ •} or all 
unary signatures. Holant (J 7 ) is tractable means that it is computable in P even when we include the description of the signatures in 
the input Q in the input size. Holant(J-") is #P-hard means that there exists a finite subset of T for which the problem is #P-hard. 
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Definition 2.4. Given a set of signatures F, we use Holant c (F) (or Pl-Holant c (F) respectively) to denote 
Holant (F U {[1, 0], [0, 1]}) ( or Pl-Holant ( F U {[1, 0], [0, 1]}) respectively). 

Replacing a signature / € F by a constant multiple c/, where c / 0, does not change the complexity of 
Holant(F). So we view / and cf as the same signature. An important property of a signature is whether it is 
degenerate. 

Definition 2.5. A signature is degenerate iff it is a tensor product of unary signatures. In particular, a symmetric 
signature in F is degenerate iff it can be expressed as X[x, y]® . 

2.2 J'-Gate and Matchgate 

A signature from J 7 is a basic function which can be used at a vertex in an input graph. Instead of a single vertex, 
we can use graph fragments to generalize this notion. An Jugate T is a tuple (H, F, tt), where H — (V,E,D) 
is a graph where the edge set consists of regular edges E and dangling edges D. Some nodes of degree 1 are 
designated as external nodes, and all other nodes are internal nodes; a dangling edge connects an internal node 
to an external node, while a regular edge connects two internal nodes. The labeling tt assigns a function from F 
to each internal node. The dangling edges define variables for the F-gate. (See Figure [l] for one example.) We 



denote the regular edges in E by 1, 2, . . . , m, and denote the dangling edges in D by m + 1, m + 2, . . . , m + n. 
Then we can define a function for this F-gate T = (H, F, tt) , 



where (yi,V2i ■ ■ ■ , Vn) S {0, 1}" denotes an assignment on the dangling edges and H(x%,X2, • ■ ■ , X m , yi, y%, . . . , y n ) 
denotes the value of the signature grid on an assignment of all edges. We will also call this function the signature 
of the F-gatc T. An F-gate can be used in a signature grid as if it is just a single node with the particular 
signature. 

Using the idea of F-gates, we can reduce one Holant problem to another. Let g be the signature of some F-gate 
r. Then Holant(FU {g}) <t Holant(F). The reduction is quite simple. Given an instance of Holant(FU {g}), 
by replacing every appearance of g by an F-gate T, we get an instance of Holant (F). Since the signature of T is 
g, the values for these two signature grids are identical. 

We note that even for a very simple signature set F, the signatures for all F-gates can be quite complicated 
and expressive. Matchgate signatures are an example. Matchgate is introduced by Valiant [551 123 |3D], whose 
definition is combinatorial in nature. Matchgates can be viewed as a special case of planar F-gates, where F 
contains Exact-One functions of all arities and weight functions ([1, 0, w], w € C) on edges. (Formally, we replace 
each matchgate edge of weight w by a path of length 2, and the new node on the path is assigned weight function 
[1,0, w].) The signature function T defined above for a matchgate is called a matchgate signature, or a standard 
signature. A signature function is realizable by a matchgate if it is the standard signature of that matchgate. 
(After a holographic transformation, a signature function is realizable under a basis if it is the transformed 
signature of a matchgate; see below.) 




Figure 1: An F-gate with 5 dangling edges. 
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2.3 Holographic Reduction 



To introduce the idea of holographic reductions, it is convenient to consider bipartite graphs. This is without loss 
of generality. For any general graph, we can make it bipartite by replacing each edge by a path of length two, 
and giving each new vertex the Equality function =2 on 2 inputs. (This is just the incident graph.) 

We use Holant (<?| 72.) to denote all counting problems, expressed as Holant problems on bipartite graphs 
H = (U, V, E), where each signature for a vertex in U or V is from Q or TZ, respectively. An input instance for 
the bipartite Holant problem is a bipartite signature grid and is denoted as 51 = (H, Q\TZ, n). Signatures in Q are 
denoted by column vectors (or contravariant tensors); signatures in TZ are denoted by row vectors (or covariant 
tensors) [T8"] . 

One can perform (contravariant and covariant) tensor transformations on the signatures. We will define 
a simple version of holographic reductions, which are invertible. They are called holographic because they 
may produce exponential cancellations in the tensor space. Suppose Holant(£|7?.) and Holant (Q'\ TZ') are two 
Holant problems defined for the same family of graphs, and T G GL2(C). We say that there is an (invertible) 
holographic reduction from Holant(C?|7?.) to Holant(C?'|7?/), and T is the basis transformation, if the contravariant 
transformation G' = T 0g G and the covariant transformation R = R'T® r map G G Q to G' G Q' and R G TZ to 
i?' G 72.', and vice versa, where G and R have arity g and r respectively. (Notice the reversal of directions when 
the transformation T® n is applied. This is the meaning of contravariance and covariance.) 

Theorem 2.6 (Valiant's Holant Theorem |40] ) . Suppose there is a holographic reduction from #G\TZ to #Q'\TZ 
mapping signature grid Q to fl' , then Holant^ = Holantfy. 

In particular, for invertible holographic reductions from Holant (<?|7?.) to Holant (Q'\ TZ'), one problem is in P 
iff the other one is, and similarly one problem is #P-hard iff the other one is also. 

In the study of Holant problems, we will commonly transfer between bipartite and non-bipartite settings. 
When this does not cause confusion, we do not distinguish signatures between column vectors (or contravariant 
tensors) and row vectors (or covariant tensors). Whenever we write a transformation as T® n F or TF, we view 
the signature or signatures as column vectors (or contravariant tensors); whenever we write a transformation as 
FT® 71 or FT, we view the signature or signatures as row vectors (or covariant tensors). 



2.4 Some Known Dichotomy Results 

In this subsection, we state some known dichotomy theorems. We first review three dichotomy theorems from 

PS- 

Theorem 2.7. Let J- be a set of symmetric signatures over C. Then Holant* (J 7 ) is computable in polynomial 
time in the following three cases. In all other cases, Holant* (J 7 ) is #P-hard. 

1. Every signature in F is of arity no more than two; 

2. There exist two constants a and b (not both zero, depending only on F), such that for every signature 
[xo, X\, . . . , x n ] G J- one of the two conditions is satisfied: (1) for every k — 0,1,..., n — 2, we have 
axk + bxk+i — axk+2 = 0; (2) n = 2 and the signature [xq,x\,X2] is of form [2aX,b\, — 2aA]. 

3. For every signature [xq, x\ , . . . , x n ] G T , one of the two conditions is satisfied: (1) For every k — 0,1, ... , n— 
2, we have Xk + = 0; (2) n — 2 and the signature [xq,xi,X2\ is of form [A,0, A]. 

The same dichotomy also holds for Pl-Holant* (J-) . 

Theorem 2.8. Let T be a set of real symmetric signatures, and letT\,T2 andTj, be three families of signatures 
defined as 

T x = {A([l,0p fc +r[0,l]® fc )|AGC,fc = l,2,...,r = 0,l,2,3}; 
T2 - {A([l,l]® fe + r[l,-l]® fe )|AGC,fc = l ! 2,...,r = 0,l,2,3}; 
J3 = {A([l,^ fe + z''[l,-^' £ )|AGC,fc = l,2,...,r = 0,l,2,3}. 

Then Holant c (J 7 ) is computable in polynomial time if (1) After removing unary signatures from J- , it falls in one 
of the three Classes of Theorem 2.7 (this implies Holant* (J 7 ) is computable in polynomial time) or (2) (Without 
removing any unary signature) J- C T\ U Ti U J-3 . Otherwise, Holant c (J r ) is #P-hard. 
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Here we explicitly list all the real signatures in J-\ U T% U J- 3, up to an arbitrary scalar factor: 



1. 


to) 


[1,0,0,..., 0,1 (or -1)], 




2. 


to) 


[1,0,1,0,...,0 (or 1)], 




3. 


to) 


[0,1,0,1,...,0 (or 1)], 




4. 


to) 


[1,0,-1,0, 1,0,-1,0, ...,0 (or 1 oi- 


-1)], 


5. 


to) 


IO, 1,0,-1,0, 1,0,-1,..., (or 1 or 


-1)], 


G. 


to) 


[1,1,-1,-1,1,1,-1,-1,... ,1 (or - 


-1)], 


7. 


to) 


[1,-1,-1,1,1,-1,-1,1,...,! (or - 


-!)]• 



Definition 2.9. A k-ary function f{x\, . . . , x k ) is affine if it has the form 

where X = (x±, X2, ■ ■ ■ ,Xk, 1), and \ is a 0-1 indicator function such that X[AX=o] is 1 iff AX — 0. Note that the 
inner product (a,X) is calculated over ¥2, while the summation over j on the exponent of i = y/—T is over ¥4. 
We use A to denote the set of all affine functions. 

We use V to denote the set of functions which can be expressed as a product of unary functions, binary equality 
functions ([1,0,1] on some two variables) and binary dis equality functions ([0,1,0] on some two variables). 

Theorem 2.10. Suppose J- is a set of functions mapping Boolean inputs to complex numbers. If J 7 C A or 
J- C V , then #CSP(J-) is computable in polynomial time. Otherwise, #CSP(J-) is #P-hard. 

As we mentioned in [TH], the class A is a natural generalization of the symmetric signatures family .F1U.F2U.F3. 
It is easy to show that the set of symmetric signatures in A is exactly J-\ U T% U J- 3. 

The following dichotomy for 2-3 regular graphs is from |27j . 

Theorem 2.11. ([27]) ^ e problem Holant([yo, y\, 2/2] |[1, 0, 0, 1]) is #F '-hard for all yo, 2/1,2/2 € C except in the 
following cases, for which the problem is in P: (1) y\ = 2/02/2; (2) j/q 2 = y\ 2 and j/02/2 = — 2/i ( Vi 7^ 0) / (3) 
yi = 0; (4) 2/0 — 2/2 — 0. If we restrict the input to planar graphs, then these four categories are tractable in P, 
as well as a fifth category j/g =2/2, an d problem remains #P-hard in all other cases. 

2.5 Characterization of Realizable Signatures by Matchgates 

A matchgate is called even (respectively odd) if it has an even (respectively odd) number of vertices. The following 
two lemmas are from [8]. 

Lemma 2.12. A symmetric signature [zq, . . . , z m ] is the standard signature of some even matchgate iff for all 
odd i, Zi — 0, and there exist r\ and r2 not both zero, such that for every even 2 < k < m, 



r\z k ~2 = r 2 z k . 

Lemma 2.13. A symmetric signature [zq, . . . , z m ] is the standard signature of some odd matchgate iff for all 
even i, Z{ — 0, and there exist r\ and r2 not both zero, such that for every odd 3 < k < m, 

rxZk-2 = r 2 z k . 

In |13) . we characterized all symmetric signatures realizable by matchgat es u nder a given basis. Here we state 
the theorem for a particular basis 



1 



1 



, which will be used in Theorem 



5.1 



Theorem 2.14. A symmetric signature [x$,x\, 
the following forms: 



, x n ] is realizable under the basis 



1 



iff it takes one of 



G 



• Form 1: there exist constants A, s, t and e where e = ±1, such that for all i, < i < n, 

Xl = X[(s + t) n - l (s - t) 1 + e(s - t) n_< (« + t)*]. 

• Form 2: there exist a constant A, such that for all i, < i < n, 

Xl = \[(n-i)(-iy + i(-iy- 1 }. 

• Form 3: there exist a constant A, such that for all i,0 < i < n. 

Xi = A[(n — 2)i]. 

3 Polynomial Interpolation 

In this section, we discuss the interpolation method we will use in this paper. Polynomial interpolation is a 
powerful tool in the study of counting problems initiated by Valiant [36 and further developed by Vadhan, Dyer 
and Greenhill 34. 20 and others. The method we use here is essentially the same as Vadhan [34] . 

For some set of signatures J-, suppose we want to show that for all unary signatures / = [x,y], we have 
Holant(J-" U {[x, y]}) <t Holant (.F). Let = (G, T U {[x, y}}, tt). We want to compute Holanto in polynomial 
time using an oracle for Holant (F). 

Let Vf be the subset of vertices in G assigned / in il. Suppose \Vf\ — n. We can classify all 0-1 assignments 
a in the Holant sum according to how many vertices in Vf whose incident edge is assigned a or a 1. Then the 
Holant value can be expressed as 

Holant = C ^V~\ (1) 

0<i<n 

where Cj is the sum over all edge assignments a, of products of evaluations at all v £ V(G) — Vf, where a is such 
that exactly i vertices in Vf have their incident edges assigned (and n — i have their incident edges assigned 1.) 
If we can evaluate these Cj, we can evaluate Holant^. 

Now suppose {G s } is a sequence of .F-gates, and each G s has one dangling edge. Denote the signature of G s 
by f s = [#sj2/s]j f° r s = 0, 1, . . .. If we replace each occurrence of / by f s in il we get a new signature grid fl s , 
which is an instance of Holant (J 7 ), with 

Holantn s = £ c^y^ 1 - (2) 

0<i<n 

One can evaluate Holantn s by oracle access to Holant (J 7 ). Note that the same set of values c, occurs. We can 
treat Ci in ^ as a set of unknowns in a linear system. The idea of interpolation is to find a suitable sequence 
{f s } such that the evaluation of Holantn, gives a linear system ([2| of full rank, from which we can solve all c,-. 

In this paper, the sequence {G s } will be constructed recursively using suitable gadgetry. There are two gadgets 
in a recursive construction: one gadget has arity 1, giving the initial signature g — [xo,yo\; the other has arity 2, 
giving the recursive iteration. It is more convenient to use a 2 x 2 matrix A to denote it. So we can recursively 
connect them as in Figure [2] and get {G s }. 







r > 




r 
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A 
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Figure 2: Recursive construction. 



The signatures of {G s } have the following relation, 
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where A 



an 012 
021 022 



and g 



We call this gadget pair (A, g) a recursive construction. It follows from Lemma 6.1 in [33] that 
Lemma 3.1. Let a, j3 be the two eigenvalues of A. If the following three conditions are satisfied 

1. det(A) ^ 0; 

2. g is not a column eigenvector of A (nor the zero vector); 

3. a/ ft is not a root of unity; 

then the recursive construction (A, g) can be used to interpolate all unary signatures. 

A similar interpolation method also works for signatures with larger arity but have two degrees of freedom. 



For example, all signatures of form [0, x, 0, y\. This is used in the proof of Lemma 4.9 



4 Dichotomy for Planar Holant c Problems 

Before presenting the main dichotomy theorem for planar Holant c problems, we prove the following theorem, 
which plays a crucial role in the proof of the main theorem. 

Theorem 4.1. Let a, b e M. 

• Ifab^l then PI- Holant? ([a, 0, 1,0,6]) is #P-hard. 

• If ab — 1 then Pl-Holant c ([a,0, 1,0,6]) is solvable in P. 

We first prove three lemmas which will be used in the proof of this theorem. 

Lemma 4.2. Let a, b,x G M, ab ^ and x ^ ±1. Then Pl-Holanf ({[a, 0,0, 0,6], [0,1,0, x]}) is #P -hard. 

Proof. Firstly, we show how to realize (=6) = [1,0,0,0,0,0,1] by [a, 0,0, 0,6]. [a, 0,0,0, 6] can be attached to a 
vertex of degree 4. We can connect 3 pairs of edges of two copies of [a, 0, 0, 0, 6] to realize the binary function 

[ fl2 ' ' 62 ]- 

If a 2 = b 2 , then we connect one pair of edges from two copies of [a, 0, 0, 0, 6] to get [a 2 , 0, 0, 0, 0, 0, b 2 }. This is 
the same as (= 6 ) = [1, 0, 0, 0, 0, 0, 1] after factoring out the non-zero factor a 2 = b 2 . 

If a 2 =/: 6 2 , then we connect [a, 0, 0, 0, 6] with a chain of [a 2 , 0, 6 2 ] of length i to get [a 2l+1 , 0, 0, 0, 6 2l+1 ]. Because 
for any i ^ j, a 2l+1 /b 2t+1 ^ a 2j+1 /6 2j+1 , we can realize (=4) — [1,0,0,0,1] using polynomial interpolation, as 
follows. Consider any signature grid on a planar graph G with n occurrences of =4 together with some other 
signatures. Let Xk/ be the sum, over all 0-1 edge assignments <r, of the products of all other vertex function 
values in G except at n vertices with =4, where k, t > and k + £ = n, and in a exactly k occurrences of = 4 have 
input 0, and exactly I occurrences of —4 have input 1. The Holant value is Ylk+£=n x k,e- Now substitute each 
occurrence of =4 by [a 2i+1 , 0, 0, 0, 6 2z+1 ]. The new signature grid has Holant value J2k+t=n x k,i(a k b e ) 2l+1 . This 
gives a Vandermonde system from which we solve for Xk t i- Now we have =4. Then we connect two copies of =4 
on one pair of edges to get =6 . 

Take a vertex of degree 6 in a planar graph attached with where the 6 incident edges are its variables. 
We will bundle two adjacent variables to form 3 bundles of 2 edges each. Then if the inputs are restricted to 
{(0,0), (1,1)} on each bundle, then the function takes value 1 on ((0, 0), (0, 0), (0, 0)) and ((1, 1), (1, 1), (1, 1)), and 
takes value elsewhere. Thus if we restrict the domain to {(0, 0), (1, 1)}, it is the ternary Equality function =3. 

Let F = [0, 1, 0, x] and let H(x\, X2, 2/1,2/2) = S 2 =n 1 F{ x iiDii Z )F{ X 2,V2, z). This H is realizable by connect- 
ing one pair of edges of two copies of F. (See Figure(3|) We will consider H as a function in (xi, £2) and (2/1,2/2)- 
However we will only connect H externally by connecting (xi,X2) and (2/1,2/2) to some bundle of two adjacent 
edges of some =6- Since = 6 enforces the values on the bundle to be either (0,0) or (1,1), we will only be interested 
in the restriction of H to the domain {(0,0), (1, 1)}. On this domain, H is a symmetric function of arity 2, and 
can be denoted as [1, 1, x 2 ]. (Note that H is not a symmetric function of arity 4 on {0, 1}, as H (0, 1, 0, 1) = x.) 

Now we have reduced Pl-Holant c ({[l, 0, 0, 1], [1, 1, x 2 }}) to Pl-Holant c ({[a, 0, 0, 0, 6], [0,1,0,2:]}). 
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Figure 3: The gadget for function H and [1, 1,2 ] 



Using (=3) = [1,0,0, 1], we can realize the Equality function =k of any arity k > 3. Then we can realize 
[1, l,x 2k ], for all k > 1. (See Figure[|) If x = 0, then we already have [1,1,0]. 

Suppose x 7^ 0. Because x 2 ^ 1 and being a positive real number, we can realize [1,1,0] by interpolation. Now 
we have reduced the problem Pl-Holant([l, 0, 0, 1] | [1,1,0]) to Pl-Holant c ({[l, 0, 0, 1], [1, 1, a; 2 ]}). The bipartite 
problem Pl-Holant([l, 0, 0, 1] | [1,1,0]) is #P-hard since it is counting Vertex Covers on planar 3-regular 
graphs g2]. □ 



The following lemma handles a special case of Theorem |4.1| The proof uses Lemma |4.2 
Lemma 4.3. Pl-HolanF([0, 0, 1, 0, 0]) is #P-hard. 



to 



Proof. We construct a reduction from Pl-Holant c ([l, 0, 0, 0, 1], [0, 1, 0, 0]), which is #P-hard by Lemma 4.2 
Pl-Holant c ([0, 0, 1, 0, 0]) by polynomial interpolation. 

Let F — [0,0, 1,0,0]. There is a series of planar gadgets (a chain of F) realizing the following sequence of 
functions: 

H 2 (x 1 ,x 2 ,y 1 ,y 2 ) = 2J F(x 1 ,x 2 ,x 3 ,X4)F(y 1 ,y 2 ,x 3 ,X4), 
£3,3:4— 0,1 

and for i > 1, 

H 2i+2 (x 1 ,X 2 ,yi,y 2 ) = H2i(x 1) X 2 ,X3,X4)H2(yi,y2,X 3 ,X4). 
£3, £4— 0,1 

The gadget for i?2i is composed of 2i functions F. As an example, the gadget for H4 is shown in Figure |4j 




Figure 4: The gadget for H4 



By calculation, 1^(0,0,0,0) = H 2i (l, 1, 1, 1) = 1, and # 2i (0, 1, 0, 1) = i? 2 i(0, 1,1,0) = H 2l (l, 0, 0, 1) = 
H 2i (l, 0, 1, 0) = 2 2l_1 , and i?2i is zero on other inputs. Again we will consider the inputs to H 2i as bundled into 
(11,2:2) and (2/1,2/2)- 

Given a planar graph G as an instance of Pl-Holant c ([l, 0, 0, 0, 1], [0, 1, 0, 0]), suppose there are n vertices in 
G attached with the function (=4) = [1, 0, 0, 0, 1]. For i = 1, 2, . . . , n + 1, we construct an instance Gi of 
Pl-Holant c ([0, 0, 1, 0, 0]) as follows: Replace each occurrence of —4 by a copy of H 2i , and replace each occurrence 
of [0, 1, 0, 0] by [0, 0, 1, 0, 0] connected with a [0, 1], which exactly realizes [0, 1, 0, 0]. Note that by replacing =4 
with H 2 i, we have bundled two adjacent edges together (in the planar embedding) for each vertex attached with 

=4- 

Let x a ,b denote the summation, over all 0-1 edge assignments a, of the products of all other vertex function 
values in G except at those n vertices with =4, where a, b > and a + b = n, and in a exactly a occurrences of 
—4 have inputs {0000, 1111}, and exactly b occurrences of =4 have inputs {0101, 0110, 1001, 1010}. 

Note that the Holant value on Gi is 

]T x ab l a (2 2l - 1 ) b . 

a-\-b—n 



On the other hand, the value of Pl-Holant c ([l, 0, 0, 0, 1], [0, 1, 0, 0]) on G is exactly x n ^. 

When we take 1 < i < n + 1, we get a system of linear equations in x a b, whose coefficient matrix is a full 
ranked Vandermonde matrix. Solving this Vandcrmonde system we obtain the value x n fi. 

□ 

The following result can be proved by interpolation as well. 

Lemma 4.4. Let a ^ {—1, 0, 1} be a real number. Then we can interpolate all [a;, 0, y, 0] and [0, y, 0, x] for x, y G C 
starting from either [0, 1, 0, a] or [a, 0, 1, 0] . 

Proof. The recursive construction is depicted by Figure [5j By a simple parity argument, every .F-gate iVj has a 




Figure 5: The recursive construction. The signature of every vertex in the gadget is [0, 1,0, a]. 



signature of the form [0, Xi, 0, yt]. After some calculation, we see that they satisfy the following recursive relation: 



Xi+l 




'3(a 2 - 


hi) a 3 + a" 






m+i_ 




3(a 3 - 


ha) a 6 + 1 




Vi 



The signatures we want to interpolate are of arity 3. But since all of them take the form [0, Xi, 0, yi] with two 
degrees of freedom, we can use the interpolation method in Section [3] Now we verify that the conditions of 

that theorem are satisfied. Let A = j?( a ., ~|~ ^ ^ a R ~!~ , , then (A, [1, a| T ) forms a recursive construction. Since 

3{a° + a) a° + 1 J 

det(A) = 3(o 4 -l) 2 7^ 0, the first condition holds. Its characteristic equation is X 2 -(a 6 +3a 2 +4)AT+3(a 4 -l) 2 = 0. 
For this quadratic equation, the discriminant A = (a 6 — 3a 2 — 2) 2 + 12(a + a 3 ) 2 > 0. So A has two distinct real 
eigenvalues. The sum of the two eigenvalues is trA = a 6 + 3a 2 + 4 > 0. So they are not opposite to each other. 
Therefore, the ratio of these two eigenvalues is not a root of unity and the third condition holds. Consider the 

second condition: if the initial vector [1, o] T is a column eigenvector of A, then we have A = A ^ , where 



A is an eigenvalue of A. From this, we will conclude that a(a 2 — l)(a 4 — 1) = 0, which can not happen given 
a ^ { — 1,0, 1}. To sum up, this recursive relation satisfies all three conditions of Lemma 3.1 and can be used to 
interpolate all signatures of the form [0, x, 0, y\. This completes the proof. □ 



Proof of Theorem 4.1 If ab = 1, then [a, 0,1,0, b) is realizable by some matchgate, by Lemma 2.12| This 
realizability also applies to the unary functions [1,0] and [0, 1]. Hence the problem Pl-Holant c ([a, 0, 1,0,6]) can 
be solved in polynomial time by matchgate computation via the FKT method 33 , 2J1 US] . In the following we 
assume that ab ^ 1 and prove that the problem is #P-hard. The case a = b = is proved in Lemma 4.3 Now 
we can assume at least one of a and b is non-zero, and by symmetry we assume a =/= 0. 

We know from our dichotomy for Holant c problems [TB] that Holant c ([a, 0, 1, 0, b]) for general graphs is #P- 
hard unless a — b=l or a = b = — 1, in which cases it is tractable. Both of these tractable cases are also included 
in the tractable cases (ab = 1) here. Therefore, if we can realize a cross function X with a planar gadget when 
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ab 1, we can reduce Holant c ([a, 0, 1, 0, b]) for general graphs to Pl-Holant c ([a, 0, 1, 0, b]) and finish the proof. 
Here a cross function X has 4 input bits, and satisfies ATqooo = ^oioi = -^loio = ^nii = 1 an d X a — for all 
other inputs a £ {0, l} 4 . 

If {a, b} <f_ {—1,0,1}, we can use Lemma [4~4] to interpolate all [x, 0,y,0], for x,y £ C. If {a, b} c { — 1,0,1}, 
then there are only four cases: [1,0,1,0,-1], [1,0,1,0,0], [—1,0,1,0,1] and [—1,0,1,0,0]. In all four cases, it is 
easy to verify that we can realize a signature with a form [ci,0,C2,0] where C1C2 ^ and c\ ^ ±C2 using the 
gadget in Figure [6j ( For [1,0,1,0,-1], we get [8,0,4,0] by using [1,0] in the gadget; for [1,0,1,0,0], we get 
[8,0,5,0] by using [1,0]; for [-1,0,1,0,1], we get [0,4,0,8] by using [0,1]; and for [-1,0,1,0,0], we get [0,1,0,3] 
by using [0,1].) After factoring out a nonzero factor, we have [c' , 0,1,0], where d £ K and c' £" {0,±1}. As a 
result, we can also interpolate all [cc, 0, y, 0], where x, y £ C. 




Figure 6: The signature of the degree 1 vertex in the gadget is [1, 0] or [0, 1]. 



Now we can use all signatures of the form [x, 0, y,0], for arbitrary x,y £ C, to build new gadgets. We also 
have all [x, 0, y] by connecting [x, 0, y, 0] to a [1, 0]. By connecting a \{ft[a, 0, \faji] to each edge of the signature 
[a, 0, 1, 0, b], we get [i, 0, 1, 0, |] for all t ^ 0, where c = ab ^ 1. Using all these, we will build a planar gadget 
in Figure [7] to realize the cross function X. In the equations below x, y,t are three variables we can set to any 
complex numbers, with t ^ 0. The parameter c is given and not equal to 1. 

(Of course we presumably could not build a cross function X if c = 1; this is exactly when the problem is in P, 
and this is also exactly when our construction of X fails. If a cross function X were to exist when c = 1 then P = 
#P would follow. However, it is still rather mysterious that algebraically c = 1 is exactly when our construction 
fails. This failure condition is by no means obvious from the equations below.) 



Figure 7: This gadget is to realize the Cross function. The signature for the center vertex (black and square) is 
[t, 0, 1, 0, ?]. The signature for the vertexes in the four corners (red and circle) is [x, 0, 1, 0]. The signature for the 
vertexes in the middle of the boundaries (green and triangle) is [y, 0, 1, 0]. 
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We can compute the signature of the gadget in Fig. [7] If the input has an odd number of Is, the value is 0. 
For other inputs, we have 

2cx 2 





Xqoqo 


4 4j. 

= x y t 


f t + 4x 3 y 2 -+ 


-4x + 


4x 2 y - 






= 2y 2 t + 


12, + y 








= -^1010 


= 2xy 2 t 


-1- 4a; 2 ?/ 2 + 4 H 


- 4xy 


2cx 
+ ^ 


Xqqh — -*1001 = -XllOO 


= -^OllO 


= x 2 y 3 t ■ 


f yt + 3x 2 y 2 


+ 3 + 


6xy + 



2 ex 

~T' 

Here we prove that for any c ^ 1, we can assign suitable complex values to x, y and t, where t 7^ 0, such that 
i = B = C ^ 0, and D — 0, where A, £?, C and D denote respectively the four functions of x, y and t listed in 
the four lines above. 

Claim 1. For any c 7^ 1, 



(x-l) 2 = ^- 
c — 1 

has a solution x tfL {0, +1, —1}. This x satisfies 



x(x + 3)\ fx + 3\ 
2 r 1 r + ex + 6 = 0. (4) 



a; — 1 / \ x — 1 

Proof. Clearly x = 1 is not a solution to (x — l) 2 = -^p Also the equation has two distinct roots. When c = 17 
there is a solution x = 2 £ {0, +1, —1}. When c 7^ 17, we can verify x — is not a solution. Hence the equation 
always has a solution other than 0, ±1. 
To verify Q we have 

(2x - 2 - x 2 - 3x)(x + 3) + (cx + 6)(x 2 - 2x+ 1) 
= -(x 3 + 4x 2 + 5x + 6) + cx 3 + (6 - 2c)x 2 + (-12 + c)x + 6 
= (c- l)x 3 - 2(c- l)x 2 + (c- 17)x 
= (c- l)x[(x - l) 2 - 16/(c- 1)] 
= 0. 

□ 

Now we fix x ^ {0, +1, —1} satisfying Q for any given c ^ 1. 
Claim 2. for any c ^ 1, we can picfc z ^ ±1 suc/i £/ia£ 

4z x(x + 3) 



(1 + z) 2 x-1 



(5) 



Proof. We are given x 7^ 0,±1. If x = —3, we can pick z = 0. Now suppose x 7^ —3. Consider the quadratic 
equation in z 

4z(x-l) =x(x + 3)(l + z) 2 . 

This is quadratic since x(x + 3) 7^ 0. We can check that z = +1 (and —1 respectively) is not a solution, as this 
would force x = — 1 (and +1 respectively). However, any solution where z ^ -1 and x 7^ 1 is equivalent to (j5b. 
Hence we have a solution z 7^ ±1 to 

□ 

Now we further fix a z 7^ ±1 satisfying ([5]), and let y — z/x such that xy 7^ ±1, for any c 7^ 1. 
Claim 3. For any c 7^ 1, there exist x (jL {0, +1, —1} and y such that xy 7^ ±1 satisfying 

2(1 + x 2 y 2 ) x + 'i 

— • + CX + 6 = 0. (6) 

(1 + xy) 2 x-1 w 
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Proof. 



2(l + x 2 y 2 ) x + 3 

—- rrr- • 7 + CX + 

(1 + xy) z X — 1 

2z \ x + 3 

2 1- ^ + cx + 6 

(1 + z) 2 / x-1 

x(x + 3)\ x + 3 
2 — I r + cx + 6 



a; — 1 7 x — 1 



= 0. 



Here we used §5§ and Q. 



Now we will set t = 4/(1 + xy) 2 . Clearly t ^ 0. We next verify that D = 0. 
By ([5) and ^ we get 

8y(l + x 2 y 2 ) 



Then 
Thus 



i ca; + 6 = 0. 
(1 + xy) 4 

t 2 y(l + x 2 y 2 ) + 2cx + 3t(l + xy) 2 = 0. 

2cx 

D = yt(l + x 2 y 2 ) + 3(1 + xy) 2 + — = 0. 



Next we show that C* = ^ 0. 



By Z? = 0, we have 



C = 2xy 2 — + 4(1 + xy) 2 - 4xy + [-yt(l + x 2 y 2 ) - 3(1 + xy) 2 ]. 

(1 + xy) 2 



Hence 



C 



8xy 2 , x2 „ 4(l + x 2 y 2 ) 

- + 1 + xy 2 - 4xy - y 

(1 + xy) 2 
2 [2xy - 1 - x 2 y 2 ] + (1 - xy) 2 

+ l) (1 - *2/) 2 
^0, 



(1 + xy) 2 

4y 
(1 + xy) 
-4y 



(1 + xy) 2 
4(1 -xy) 2 



1 - x 



using 

The next task is to show B = C. 
We have 



Hence 



B 



4(1 - xy) 2 
1 - x 



4(1 -xy) 



C = 4(1 - xy) 2 + xS 
4(1 -xy) 2 



1 -x 



4(1 - xy) 5 
1 - x 



C. 



Finally we verify A = C as well. 



□ 



A = (x 4 y 4 + 1)* + a[C - 2xy 2 t] = C+ (x - 1)C - 2x 2 y 2 i + (x 4 y 4 + l)t = C - 4(1 - xy) 2 + t(x 2 y 2 - l) 2 = C. 



□ 



Now we come to the main dichotomy theorem for Pl-Holant c problems. 
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Theorem 4.5. Let J- be a set of real symmetric signatures. Pl-Holant c (J r ) is #P-hard unless J- satisfies one of 
the following conditions, in which case it is tractable: 

1. Holant c (J r ) is tractable (for which we have an effective dichotomy or 

2. Every signature in T is realizable by some matchgate (for which we have a complete characterization JEj). 

Before we give the proof, we do some normalization of the signature set J 7 . Since any degenerate signature 
[a;,y]® fe can be replaced by the corresponding unary signature [x,y] without changing the complexity of the 
problem, we always assume that all the signatures in J 7 , whose arity is greater than 1, are non-degenerate. Since 
[1,0] and [0, 1] are freely available, we can construct any sub-signature of an original signatures as well as any 
signature realizable by some .F-gate. 

The main idea of the proof is to interpolate all unary functions. If we can do that, we can reduce the problem 
Pl-Holant* ( F) to Pl-Holant c (J-") and finish the proof. We note that our dichotomy in [IB] for Holant*(J r ) also 
holds for planar graphs. In some cases, we cannot interpolate all unary functions, then we prove the theorem 
separately, mainly using Lemma 4.2 and Theorem |4.1| The following lemma is for interpolation of unary functions. 



Lemma 4.6. If we can construct from J- a gadget with signature [a, b,c], where b 2 ^ ac, b ^ and a + c ^ 0, 



then we can interpolate all unary functions. (Hence the conclusions of Theorem 4-5 hold.) 



Proof, we use the interpolation method as described in Section [3] We consider two recursive constructions 



a b 




T 


) and ( 


a b 




¥ 


b c 


1 





b c 




l 



We use A to denote 



), and argue that at least one of them will succeed given the conditions on a, b, c. 
Since b 2 ^ ac, A is non-degenerate, the first condition of Lemma 3.1 



is satisfied for 



both recursive constructions. If both [1, 0] and [0, 1] are column eigenvectors of A, then b = 0, a contradiction. 
So at least for one of the two recursive constructions, the second condition of Lemma [3. 1| is satisfied. Since A is 
a real symmetric matrix, both its eigenvalues arc real. If the ratio of two real numbers is a root of unity, they 
must be the same or opposite to each other. If the two eigenvalues are the same, we have 6 = and a — c, a 
contradiction. If the two eigenvalues are opposite to each other, then we have a + c = 0, also a contradiction. 



Therefore, the third condition of Lemma 3.1 is also satisfied for both recursive constructions. To sum up, at least 
one of the two recursive constructions satisfies all the conditions of Lemma |3.1[ This completes the proof. □ 



If we can construct from J- a gadget with a binary symmetric signature [a, 6, c], which satisfies all the conditions 
in Lemma |4.6| then we are done. For most cases, we prove the theorem by interpolating all unary signatures. 
However, in some more delicate cases, we are not able to do that. For example, if all signatures from T have 
the parity condition, which includes a proper superset of matchgate signatures, then all unary signatures we can 
realize have form [a, 0] or [0, a], so we can not interpolate all unary signatures. For these cases, our starting point 
is Theorem 14.11 

We define some families of symmetric signatures, which will be used in our proof. 
Gi = {[a, 0,0,..- ,0,6] | a6^0} 

G2 = {[xq, %i, ■ ■ ■ , %k] I Vi is even, = or Vi is odd, = 0} 

03 = {[xo,Xi, ■ ■ ■ ,X k ] | Vi, Xi + X i+2 = 0} 

■M = { f | / is realizable by some matchgate }. 

We note that Gi , G2 and G3 are supersets of J-\ , T2 and J-3 respectively. Furthermore (the real part of) T% C 
■M. C Qi. The conditions in G2 are called parity conditions. The following several lemmas all have the form "If 
J- % A, then the conclusions of Theorem |4.5| hold." After proving each lemma, in subsequent lemmas, we only 
need to consider the case that F C A. 



Lemma 4.7. If J- % Gi U G2 U G3, then the conclusions of Theorem 



hold. 



Proof. Since T % Gi U G2 U ^3, there exists an / 6 T and /^SiU^U^. Since all unary signatures are in G3, 
the arity of / is greater than 1 and / is non-degenerate. There are two cases according to whether / has a zero 
entry or not. 
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(1) / has some zero entries. If there exists a sub-signature of / of the form [0, a, b] or [a, 6,0], where ab 7^ 0, 
then we are done by Lemma |4.6| Otherwise, we can conclude that there are no two successive non-zero entries. 
So the signature / has this form [0*°xi0 ll X20 12 • • • XfcO Ifc ], where k > 1, Xj ^ and for all 1 < j < k — 1, ij > 1. 
If for all 1 < j < k — 1, ij is odd, (including k = 1), then / G G2, a contradiction. Otherwise there exists a 
sub-signature of form [x, 0, 0, • • • , 0, y], where xy 7^ and there are a positive even number of 0s between x and 
y. If this is the entire /, then / G Gi, a contradiction. So there is one before a; or after y. By symmetry, we 
assume there is a before x, so we have a sub-signature [0,2,0,0, • • • ,0,y], whose arity is even and at least 4. 
We label its dangling edges 1, 2, ■ • • , 2k. Then for every i = 1, 2, • • • , k — 1, we connect dangling edges 2i + 1 and 
2i + 2 together to form a regular edge. After that, we have an .F-gate with arity 2, and its signature is [0, x, y\. 
Then we are done by Lemma |4~6} 

(2) / has no zero entry. We only need to prove that we can construct a function [a', b', c'] satisfying the three 
conditions in Lemma |4~6} Suppose all sub-signatures of / with arity 2 do not satisfy all the three conditions. For 
each sub-signature [a', &',</], either a' + c' = 0, or b' 2 = a'c'. If all of them satisfy a' + d = 0, then / G Q 3 . A 
contradiction. If all of them satisfy b' 2 — a'c' , then / is degenerate. A contradiction. W.l.o.g., we can assume there 
is a sub-signature [a, 6, c, d] of /, such that a + c = 0, b + d 7^ 0, and c 2 = bd. We get this sub-signature [a, 6, c, d] by 
[1, 0] and [0, 1]. Combining two [a, b, c, d], we can get a function [a', b' , c'] = [a 2 +2b 2 +c 2 , ab+2bc+cd, b 2 +2c 2 +d 2 ] — 
[2(b 2 + c 2 ),c{b + d),(b + d) 2 }. Then b' = c(b + d)^0. a' + c' > 0. And a'c' - b' 2 = {b + d) 2 {2b 2 + c 2 ) > 0. We are 
done by Lemma |4.6| □ 



The following lemma uses Theorem |4.1| in an essential way, which in turns depends on the crossover. 
Lemma 4.8. // J- % Gi U M. U G3, then the conclusions of Theorem \4-.5\ hold. 



Proof. If T % <?iU£ 2 U<?3, then by Lemma [4~7| we are done. Otherwise, there exists a signature / G T C £iU£ 2 U<?3 
and / ^ Gi U M U £3. Then it must be the case that / € G%. Note that every signature with arity at most 3 in G2 
(this is called the parity condition) is also contained in M, so / is of arity greater than 3. Let / = [xq, xi, ■ ■ ■ , x n ], 
for some n > 4. Suppose there exists some i & [2, 3, • • • , n — 2] such that Xi 7^ 0. If Xi^ 2 %i+2 ^ x h then we can 
get \xj-2, 0, Xj, 0, Xi+2 \ by [1,0] and [0,1] which restrict the signature to a sub-signature. Then the problem is 



f^P-hard by Theorem 4.1 and we are done. Otherwise, we have Xi- 2 Xi+ 2 — xf 7^ 0. Then starting from Xi- 2 7^ 
and if z — 2 G [2, 3, • • • , n — 2], we can get x^^Xi — x 2 _ 2 7^ 0. Similarly we can start with x i+2 . A signature 
satisfying the parity condition and is a geometric series on the alternate entries is realizable by a matchgate [71 [S] , 
a contradiction. 

Now we may assume Xi = for all i € [2, 3, • • • , n — 2]. Since f E G2 — (-M U Gi), we know that there are 
only three possible subcases: (1) n is odd, n > 5, 7^ and si = i n = 0; (2) n is odd, n > 5, 7^ 

and = 2>n-i = 0; (3) n > 6 is even, xia;„_ 1 7^ and xq — x n = 0. This uses the theory of matchgate 
realizability [7J|S]. Crucially, if n is even and n < 6, then 71 = 4 and the case xix n ~i 0, xq = x n = belongs 
to M.. The subcases (1) and (2) are reversals of each other and (3) contains a signature in form (1) and (2). So 
after normalizing (and connecting pairs of edges together if n > 5), we will get a signature [0, 1, 0, 0, 0, x] where 



x 7^ 0. So we have both sub-signature [0, 1,0,0] and [1,0,0,0,2;]. As we proved in Lemma 4.2 the problem is 



#P-hard and we are done. This finishes the proof. □ 



Lemma 4.9. If [0,1,0, a;] G J- (or [l,0,a;,0] G J-) where x G M., x 7^ ±1, then the conclusions of Theorem 4-5 
hold. 



Proof. If x 7^ 0, we can use Lemma 4.4 to interpolate [0,1,0,0]. So we assume we have [0,1,0,0] from T . If 
J~ % Gi U M. U G3, then by Lemma |4.8| we are done. If J- C A4, then the problem is tractable and we are done. 
Otherwise, there exists a signature f G T C Gil) MU G3 and f ^ M. That is / G (Gi U £3 - M). 

If / has arity > 1 and of the form [xq, xi, —xq, —x\, Xq • • •] G G3 — then we will have xqX\ 7^ 0. Otherwise 
we would have / G M, a contradiction. Connecting one unary signature [xo,xi] to [0,1,0,0], we get [xi,xo,0] 
which satisfies all the conditions in Lemma |4.6[ and we are done. 

Now we consider / = [1, 0, 0, • • • , 0, y] G Gi — M., where y 7^ 0. Since / g" /A, its arity n is greater than 2. If n 
is odd, we can connect its edges except one to get a unary signature [1, y]. Then we can use a similar argument as 
above and we are done. If n is even, then it is at least 4, since / S - M.. After connecting its edges except four, we 



can get [1,0, 0,0, y]. Together with [0, 1,0,0], we know the problem is #P-hard by Lemma 4.2 This completes 



the proof. □ 
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Lemma 4.10. If J- <2 Gi U T% U G3, then the conclusions of Theorem\4-5\ hold. 



Proof. If J- ' % ^iU^2U^/3, then by Lemma 4.7 we are done. Otherwise, there exists a signature / £ T C Q1UQ2UQ3 
and / ^ 5i U J2 U 63. Then it must be the case that f € Gs- Note that every signature with arity less than 3 in 
t/ 2 is also contained in Gi U ^3, so / is of arity greater than 2. Since f Gi, there is some non-zero in the middle 
of the signature /, after normalization, we can assume there is a sub-signature of form [0, 1,0, x] (or [x, 0,1,0]). 
If x 7^ ±1, then by Lemma |4.9| we are done. Otherwise, for every such pattern, we have x = ±1. Since / £ 
then there is some sub-signature [0, 1, 0, —1] and because f $ G3, there is some sub-signature [0, 1, 0, 1]. Therefore, 
there is a sub-signature [1, 0, 1, 0, —1] of / 
are done. This completes the proof. 



Then by Theorem 4.1 we know that the problem is #P-hard and we 

□ 



Lemma 4.11. // J- $2 Gi U G3, then the conclusions of Theorem 4-5 hold 



Proof. If T % Gi U J-2 U ^3, then by Lemma [4.10| we are done. Otherwise, there exists a signature / € T C 
Gi U T2 U G3 and / £ Gi U £3. Then it must be the case that / £ J r 2. Note that every signature with arity less 
than 3 in is also contained in Gi U £3, so / is of arity at least 3. Then / has a sub-signature [1,0,1,0] or 
[0, 1, 0, 1]. By symmetry, we assume it is [1, 0, 1,0]. If T C T\ U T2 U J3, then Theorem 4.5 trivially holds and 
there is nothing to prove. If not, there exists a signature j e J - Ji U J 2 U J3. By T C U .F2 U £3, either 
.g G Si - 7i U T 2 U J" 3 (C Gi - Fi) or 5 e g 3 - ^"i U J- 2 U J- 3 (C G3 - H). 

For the first case, g £ (Gi — after a scale, g is of form [1, 0, 0, • • • ,b], where b ^ {— 1, 0, 1}. If the arity 
of g is odd, we can realize [1,6]. (We connect every two adjacent dangling edges into one edge and leave one 
dangling edge.) Then connecting this unary signature to one dangling edge of [1,0, 1,0], we can realize a binary 



signature [1,6, 1]. Then by Lemma 4.6 Theorem 4.5 holds. If the arity of g is even, we can realize [1,0,6] (leave 
two dangling edges). By connecting one of its dangling edge to one dangling edge of [1, 0, 1, 0], we can have a new 
ternary signature [1, 0, 6, 0]. By Lemma 4.9 we are done. 

For the second case g £ (G3 — ^3), g has a sub-signature of form [1,6], where 6 ^ { — 1,0,1}. By the same 
argument as above, Theorem |4.5| holds. This completes the proof. □ 



Lemma 4.12. If ' T % Gi 1 ^ T3, then the conclusions of Theorem\4-5\ hold. 



Proof. If T % Gi U G3, then by Lemma |4.11| we are done. Otherwise, there exists a signature / £ T C Q 1 U G3 
and / Q\ U T3. Then it must be the case that / £ G3, and / has a sub-signature of form [l,a, — 1], where 
a £{-1,0,1}. 



If T C {[1, 0, 1]}U03, then Holant* (J-) is polynomial time computable by Theorem 2.7 and as a result Theorem 
4.5 trivially holds and we are done. 



If not, there exists a signature g £ J- C Q x U G3 and g £ {[1, 0, 1]} U Q3. Then it must be the case that g £ Gi- 
The arity of g is greater than 1, as g £ C/3. 

If the arity of g is 2, then g is of form [1, 0, 6], where 6 £ {—1, 0, 1}. Connecting two signatures [1, 0, 6] to both 
sides of one binary signature [1, a, —1], we can get a new binary signature [1, ab, — b 2 ]. It satisfies all the conditions 
of Lemma 4.6 and we are done. If the arity of g is greater than 2, then we can always realize a signature [1, 0, 0, 6], 
where 6 7^ 0. (We connect the unary signature [1, a) to all its dangling e dges except the three ones.) Then we 
can use an .F-gate in Figure [8j Its signature is [l,a 2 6, 6 2 ], and by Lemma 4.6 we are done. This completes the 




Figure 8: The function on degree 2 nodes is [1, a, —1], and the function on degree 3 nodes is [1, 0, 0, 6]. 



proof. 



□ 



By the above lemmas, the only case left we have to handle is that J- C Q 1 U J-3. This is done by the following 
lemma, which completes the proof of Theorem |4.5| 
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Lemma 4.13. If J 7 C Gi U J-3, then the conclusions of Theorem\4-5\ hold 



and then by Theorem 2.7 part (3), Holant c (J r ) is 
Now suppose J ^ Ji U J3 and 



Proof. If J 7 C T\ U J-3, then by Theorem 2.8 part (2), Holant c (J r ) is computable in polynomial time. Similarly, 
if J 7 C WUJ 3 U{[1,0, 1]}, then by Theorem [2T8] part (1), 
computable in polynomial time. Hence in these two cases, Theorem |4.5| holds. 
J2MUJ3U{[1,0,1]}. 

There exists / e J 7 — Ji U J3. Since J 7 C ^ U J3, such an / € 

Now there are two cases. The first case is that we have such an / $.U, and so, / £ Jnft- (J-\ U J3UW). 
The arity of / is greater than 1. By connecting its dangling edges together except two or three depends on the 
parity of the arity of /, we can assume / has form [1, 0, a) or [1, 0, 0, a], where a ^ {—1, 0, 1}. 

The second case is every / e J 7 n G\ — {T\ U J-3) is also in U. By J ^ WU J 3 U {[L 0, 1]}, there exists 
fx G T - (U U J3 U {[1,0, 1]}). Since J 7 C Q x U J" 3 , and /1 g J" 3 , we get /1 G Gi- If /1 £ Fi, we could use this 
/1 as the / above, namely f\ E J 7 <~) Gi — (J~i U J3 U W). A contradiction. Thus f\ G J-j.. Also we have some 
/2 G J 7 — (J-i U J-3). So / 2 e 5i, since J 7 C (/! U J^. Also since we are in this second case, certainly f 2 G W. 

So we have fx, / 2 G J 7 n £1 such that /1 G 7"i but /1 ^WU7 3 U{[1, 0, 1]}, and f 2 € U but / 2 G" J 7 !. The 
arity of /1 is at least 2. We claim it is greater than 2. Otherwise, /1 being from J 7 ! and not [1, 0, 1], it would be 
fx = [1, 0, —1] G J-3, a contradiction. So f\ has form [1, 0, 0, ... , ±1] of arity at least 3. fi is of form [1, a'], where 
a' { — 1, 0, 1}; this follows from / 2 G U n £1 — J-"i. By connecting all the dangling edges of f\ except two with 
/2, we can construct an Jugate with signature of form [1,0, a], where a ^ { — 1,0,1}. This is one of the above 
two forms after the first case. To sum up, in both cases, we have some / of the form [1, 0, a] or [1, 0, 0, a], where 
a £ {-1,0,1}. 

If T C Gi U {[0,1,0]} UW, then by Theorem [2~8| part (1), and then by Theorem [27] part (2) (with a = 
and b = 1), Holant c (J-") is computable in polynomial time and Theorem 4.5 holds. Otherwise, there exists 
g G J- C Gi U J-3, and g C?i U {[0, 1,0]} U Then 5 must be in and have one of the following sub- 
signatures: [1, 1, —1], [1, —1, —1], [1, 0, —1, 0], [0, 1, 0, —1]; this follows from a careful examination of the forms of 
J-3. By symmetry (taking the reversal of both / and g) 7 we only need to consider the cases / = [1,0, a] or 
[1,0,0, a], where a <£ {-1,0,1}, and g = [1,1,-1] or [1,0,-1,0]. 

According to / and g, we have four cases. If / = [ 1, 0, a] and g = [1, 1, —1], then connecting them together 
into a chain fgf, we can realize [1 



4.6 



we are done. If / = [1 
Then we can realize [1,0,- 



0,o] and g = [1,0,-1,0] 



, 0]. So by Lemma 



4.9 



for 

we 



, a, —a ]. By Lemma 
each dangling edge of g, we extend it by one copy of /. 
are done. If / = [1, 0, 0, a] and g — [1, 1, —1], we can connect a unary signature [1, 1] (sub-signature of g) to one 
dangling edge of /, and realize a binary signature / = [1,0, o]. This reduces it to the first case, which has been 
proved. If / = [1, 0, 0, a] and g = [1, 0, — 1, 0], we can realize a unary signature [1, a] from / by connecting two of 
its dangling edges together, and then connect this unary signature to one dangling edge of g to realize [1, —a, —1]. 
Note that [1, —a, —1] Gi U J-3, by Lemma 4.12 we are done. □ 



5 Dichotomy for Planar Weighted #CSP 



In this section, we prove a dichotomy for planar real weighted #CSP. Compared to the dichotomy for general 
real weighted #CSP, the new tractable cases for planar structures are precisely those which can be computed by 
holographic algorithms with matchgates. Since all the equality functions are assumed to be available, the only 

"1 1 " 



possible basis used in holographic algorithms is 



1 



-1 



(this can be computed by the characterization in |11|). 



Now we present the dichotomy theorem for planar weighted #CSP. 

Theorem 5.1. Let J- be a set of real symmetric functions. P1-#CSP(J 7 ) is #P-hard unless J 7 satisfies one of 
the following conditions, in which case it is tractable: 



1. #CSP(J 7 ) is tractable (for which we have an effective dichotomy J16jj): or 

1 



2. Every function in T is realizable by some matchgate under basis 
characterization JBI)- 



1 



(for which we have a complete 



The main proof idea is to reduce Pl-Holant c problems to P1-#CSP problems. P1-#CSP(J") is exactly the same 
as planar Holant with all the Equality functions, i.e., Pl-Holantp 7 U {[1, 1], [1, 0, 1], [1,0, 0, 1], [1,0,0, 0, 1], . . .}). 
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We can use a holographic reduction under the basis H 



1 1 
1 -1 



Under this transformation, the problem is 



transformed to, and hence has the same complexity as Pl-Holant(i?J r U {[1, 0], [1, 0, 1], [1, 0, 1, 0], [1, 0, 1, 0, 1], . . .}). 
Since this holographic reduction gives us [1,0] (from [1, 1]), if we can further realize (or interpolate) [0, 1], we can 



view the problem as a Pl-Holant c problem and apply Theorem 4.5 to HTU {[1, 0, 1], [1, 0, 1, 0], [1, 0, 1, 0, 1], . . .} to 



get a proof of Theorem 5.1 In the following, we show how to realize (or interpolate) [0, 1]. Once we have [0, 1], 
the translation of the criterion of Theorem |4.5| to Theorem |5.1| is straightforward. 

It turns out that to realize (or interpolate) [0, 1] in some cases is difficult. The following lemma says that 
it is also sufficient if we can realize (or interpolate) [0,0,1]. [0,0,1] can be viewed as two copies of [0,1], as 
[0,0, 1] = [0, 1] <g> [0, 1]. Intuitively, we will use one copy of [0,0, 1] to replace two occurrences of [0, 1]. However, 
there are two technical difficulties. One is that there may be an odd number of occurrences of [0, 1] used in the 
input instance; the second difficulty, which is more subtle, is that we have to pair up two copies of [0, 1] while 
maintaining planarity of the instance. 

Lemma 5.2. Pl-Holant{F U {[1, 0], [0, 0, 1], [1, 0, 1], [1,0,1,0], [1,0,1,0,1],...}) is #P-hard (or in P) if and only 
if Pl-Holant c {F U {[0, 0, 1], [1, 0, 1], [1, 0, 1, 0], [1, 0, 1, 0, 1], . . .}) is #P-hard (or in P). 

Proof. There is one more function [0, 1] in the second signature set than the first, so obviously the first one can 
be reduced to the second one. Hence if the second problem is in P, so is the first. We have already proved a 
dichotomy theorem for Pl-Holant c problems. So now we may assume the second problem is ^P-hard, and show 
that the first problem is also #P-hard. 

We observe that all the proofs in this paper and |16j , when the second problem for any signature set is proved 
to be #P-hard, one of the following three problems: (a) Pl-Holant([l, 0, 0, 1]|[1, 1,0]), (b) Pl-Holant([l, 1,0,0]), or 
(c) Holant[0, 1,0,0] (respectively counting Vertex Cover, Matching for planar 3-regular graphs, or Perfect 
Matching for general 3-regular graphs) is reduced to it by a chain of reductions. There are only three reduction 
methods in this reduction chain, direct gadget construction, polynomial interpolation, and holographic reduction. 

Given an instance G of Pl-Holant([l, 0, 0, 1]|[1, 1,0]), Pl-Holant([l, 1,0,0]), or Holant[0, 1,0,0], we consider the 
graph G U G, which denotes the disjoint union of two copies of G. 

Notice that the value of Pl-Holant([l, 0, 0, 1]| [1, 1, 0]), Pl-Holant([l, 1, 0, 0]), or Holant[0, 1, 0, 0] on the instance 
G is a non-negative integer, and the value on G U G is its square. So we can compute the value on G uniquely 
from its square. Suppose the reduction chain on the instance G produced instances G\, G%, . . . , G m of the second 
problem. The same reduction applied to GU G produces instances of the form Gi U G\, Gi U G2, . . . , G m > U G m / . 
(We note that the reduction on GU G may produce polynomially more instances than on G because of polynomial 
interpolation.) 

Now we only need to show how to transform instances Gi U G\ , G2 U G2 , . . . , G m ' U G m i in the second problem, 
to instances of the first problem with the same values (replacing all occurrences of the signature [0, 1] by some 
[0, 0, 1]). Gi U Gi is a planar graph with zero or more vertices of degree one attached with the function [0, 1]. We 
want to use one copy of [0,0, 1] to replace one pair of [0, 1], while maintaining planarity. 

Take a spanning tree of the dual graph of Gi. Let the outer face be the root. Choose an arbitrary leaf of 
this tree, which corresponds to a face G of Gi. Suppose C is the face corresponding to the parent of G in the 
tree. If there are an even number of vertices of degree one attached with [0, 1] in face G, we can perfectly match 
them and realize them using [0, 0, 1] while maintaining planarity in this face. This can be done by matching these 
dangling vertices of degree one in a clockwise fashion on this face G. If there are an odd number of [0, 1] in face 
G, we choose one edge e between G and C , and add a new vertex v e on e, and connect two new vertices of degree 
one to v e . The two new vertices are attached [0, 1], and v e has degree 4 and is attached [1, 0, 1, 0, 1]. The effect 
of [1, 0, 1, 0, 1] connected by two [0, 1] is the same as the function [1, 0, 1], which is exactly the same as the edge e 
itself. We put one new vertex with [0, 1] in face G, and the other one in face G'. Now, there are an even number 
of [0, 1] in face G, and we can replace them by [0, 0, 1] in G, as before. We may repeat this process, until we reach 
the root in the dual graph of Gi. If we do the same for the two Gi in G^ U Gi, we will have an even number of 
[0, 1] in the common outer face and can at last perfectly match the [0, 1] vertices and realize them by [0, 0, 1]. In 
the end we get an instance of the first problem, which has the same value. □ 



To sum up the above discussion, and apply Theorem 4.5 we have the following lemma, which is the starting 
point of our proof of Theorem |5.1| 



18 



Lemma 5.3. If we can realize (or interpolate) [0, 1] or [0,0, 1] from H JU{[1, 0], [1,0,1], [1,0, 1,0], [1,0, 1,0, 1],. . .}, 
then the conclusion of Theorem\5.1\ holds. 



Next we give two lemmas which give a general condition to realize or interpolate [0, 1] or [0,0, 1]. 

Lemma 5.4. Let a £ K. If a $ {0, 1, —1}, then we can interpolate [0, 1] from ({[1, a], [1, 0], [1, 0, 1], [1, 0, 1, 0], . . .}). 

Proof. For every j > 1, we can take a function Fj + i = [1, 0, 1, 0,1,.. .] of arity j + 1, and connect j functions [1, a] 
to it. The row vector form of the function (i.e., a listing of its values) of arity j composed of j copies of [1, a] is 



The column vector form of Fj + \ isl/2^^J +1/2 y ^ j . The 2 J x 2 matrix form of Fj + \ is 

1/2 ( \ )* J ®(l,l) + l/2( _\ Y J 0(1,-1). 
Our gadget realizes 



1/2 I J ) ® (1,1) + 1/2 ( _\ J 



^(M) + ^d,-i). 



Because a € K and a ^ {0, 1, —1}, (1 + a)/(l — a) is well defined and is neither zero nor a root of unity. We 
can interpolate any unary function x(l, 1) + y(l, —1), in particular [0, 1]. □ 

Lemma 5.5. Let a € K. 7/ a $ {0, 1, — 1}, then we can interpolate [0, 0, 1] from [1, 0, a] . 

Proof. The function of a chain of length j composed of [1, 0, a) is [1, 0, a- 7 ']. Since the real number a ^ {0, 1, —1}, 
we can interpolate all [x, 0, y], and in particular [0, 0, 1], by polynomial interpolation. □ 



Proof of Theorem |5.1[ In this proof, we augment the class T\ U Ti U -F 3 to include those degenerate signatures 
which can be obtained from tensor products from unary signatures in T\ U T% U J3. 

If 77 J 7 C J 7 ! U J-2 U J-3 , then the problem is tractable even for general graphs and the conclusion of the theorem 
holds. Now we assume that there exists an / G HT — (T\ U Ti U J-3). In the following, we will prove that we can 
realize (or interpolate) [0,1] or [0,0,1] from / and {[1,0], [1,0,1], [1,0,1,0], [1,0,1,0,1],...}. 

The general thrust of the proof is to squeeze all possible / into several standardized forms, and either prove 
f^P-hardness or reach a contradiction. We assume for a contradiction that we cannot realize (or interpolate) [0, 1] 
or [0, 0, 1]. Suppose / = [/o, /1, • • • , /n]- Since we have [1, 0], we can always take an initial subsequence of an / 
we already have as the signature of a realizable function. Given a function g with arity r > 1, we often use the 
gadget composed of two copies of g such that r — 1 inputs of them are connected to each other. We call this 
the double gadget from g. We separate two cases according to whether f Q = 0, or / 7^ which we normalize to 
/o = l- 

1. fo = 0. 

As the constant function is in T\ U T% U J-3, / is not identically 0, and thus for some i > 1, /j ^ 0. If 
fo = and f\ ^ 0, then we can connect n — 1 functions [1, 0] to / to get [0, /1], which is [0, 1] up to a factor. 

So we have /o = fi = 0, then n > 2. If f^ ^ 0, then we can connect n — 2 functions [1,0] to / to get 
[0, 0, /y, which is [0, 0, 1] up to a factor. 

So we have /o = /1 = fa = 0, then n > 3. Let m < n be the first nonzero, fo = /1 = f% = ■ ■ ■ = f m -i = 0, 
/,„ 7^ 0, then we can connect a function [1,0, 1] to two dangling edges of / to get a function whose first 
nonzero entry is f m at index m — 2. We can repeat this process until exactly one or two zeros are left at 
index or at index and 1, and we reach one of the two scenarios above. 

2. fo = I- 



By Lemma 5.4 we only need to consider fx g {0, 1, —1}. Otherwise, we are done. 
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(a) fo = 1 and fa = ±1. 

If n = 1, then / = [fa, /i] = [1, ±1] 6 7i U J2 U ^3, a contradiction. Therefore we have n > 2, we can 
take its initial part [1, fa, fa]. Connecting one edge to [1, fa], we get [1 + /j , fa + /1/2] = [2, ±(1 + / 2 )]. 



By Lemma 5.4 we only need to consider fa <E {1, — 1, —3}. 

We can construct another gadget which connects two inputs of [1, 0, 1, 0] by [1, /1, fa]- This produces a 
unary signature [I + /2, 2/i]. It follows that fa ^ — 1, since otherwise we have [0, 1] after normalization. 
Next we rule out fa = —3. The double gadget of [fa, fa, fa] = [1, ±1, iy has signature [2, —2, 10] and 
[2, 2, 10]. After normalizing, this gives [1, ±1, 5] and 5 {1, — 1, —3}. Hence, we may assume fa = 1. 
Our goal in this case 2. (a) of fa = 1 and /1 = ±1 is to extend this pattern [1, ±1, 1, . . .]. Assume we have 
proved that fa = 1 (respectively fa — (— l) 3 ) for j = 0, 1, . . . , m (m > 2). Since / ^ Ji U J2 U ^"3, the 
arity n > m. We can connect [/ ,/i] to [/o,/i,- • ■ ,/m+i], to get a function [fo+f?Jofa + fafaJofa + 
/1/3, . . .] of arity m, which is [2, 2, . . . , 1 + f m +i] (respectively [2, -2, . . . , - f m , f m - f m +i})- By 
what has been proved inductively, f m +\ = 1 (respectively f m +i = fm-i)- So in this case we showed 
that either / e Ji U J2 U J3, which is a contradiction, or #P-hardness. 

(b) f a = 1 and /1 = 0. 

Since [1,0] £ .F2, and f & J-2, we have n > 1. If / were degenerate it would be [1, 0]®" = [1, 0, . . . , 0], 
which would belong to the augmented class of J-\ U J- 2 U J- 2 . But / does not. So / is non-degenerate, 
in particular there is some nonzero entry other than fo. Suppose f m is the first nonzero f, for i > 1. 
Because we can connect some [1, 0, 1] to [/ , 0, . . . , 0, f m ] to get [1, f m ] or [1, 0, f m ], we have f m = ±1 
by Lemma |5 . 4| and 1 5 . 5 1 Since / ^ T\ we have n > m. 

We prove m is an even number. Otherwise, we can get [1, f m , f m +i] by connecting some [1,0, 1]. By 
the proof for case 2. (a), we get f m +\ = 1- We can also get [1, 0, 0, f m , / m +i], since rn > 1, whose 
double gadget has the signature [1 + fa, f m f m +i, 3 fa + fa+i\ = [2, ±1,4]. This gives #P-hardness, 
as we can get [2, ±1]. 

Now we know m must be even. Next we show that in fact m = 2. Otherwise, m > 4 and we 
can get [1, 0, 0, 0,f m , f m+ i], whose double gadget has the signature [1 + fa, fmfm+iAfm + fa+i\ = 
[2,±/ m+ i,4 + fa + i\- By what has been proved so far this also leads to #P-hardness. 
Now we have m — 2 and reached [1, 0, ±1, / 3 ], whose double gadget has the signature [2, ±/ 3 , 2 + /|], 
./;; 0. 

Again since / ^ Ji U J2 U J3 we have n > 3. Hence we have [1, 0, ±1, 0, /J. For [1, 0, —1, 0, /y, by 
connecting two edges with [1, 0, 1], we get [0, 0,-1 + f&], and we must have fi = 1, or else we have 
the signature [0, 0, 1], after normalization. For [1, 0, 1, 0, fa], by connecting two edges with [1, 0, 1], we 



get [2,0,1 + fa], and it follows from Lemma 5.5 that fa 6 {1,-1,-3}. Connecting three edges of 



[1, 0, 1, 0, fa] to three edges of [1, 0, 1, 0, 1], we get [4, 0, 3 + fa], which rules out fa = —1, by Lemma 5.5 
again. The double gadget of [1, 0, 1, 0, fa] gives [4, 0, 3 + /|], which rules out fa = —3. To sum up, we 
get fa = l. 

We have reached [1, 0, ±1, 0, 1, . . .]. The rest of the proof is similar to the induction proof for the case 
2. (a) but by skipping all entries with an odd index. Assume we have proved that fj are of the proper 
form, for j = 0, 1, . . . , m. More precisely, fa — for all odd j < m, and, either f 2 j = 1 for all j < m/2, 
or f 2 j = ( — iy for all j < m/2. Since / ^ T\ U T 2 U Tz, the arity n > m. We connect [fa, fi, fa] to 
[fa, fa, ... , f m+ i], to get a function g of arity m — 1. If m is even, f m —i — 0, and f m +i is added to or 
subtracted from f m -i, namely f m -i ± fm+l, to form the last entry in g at index m — 1. This entry 
should be zero by induction, so f m +i =0. If m is odd, we can repeat the proof in the case 2. (a), but 
we ignore all zero entries at odd indexed locations, then the induction can be completed as before. 
This completes the proof. 

□ 



6 Dichotomy for Planar 2-3 Regular Graphs 

In this section we prove a dichotomy for Holant on planar 2-3 regular graphs. This setting is very interesting for 
at least two reasons. From dichotomy theorem point of view, this is the simplest nontrivial setting and always 
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serves as the starting point of more general dichotomy theorems as in [16l [9] . This was also a focus of several 
previous work Q31 [551 HH1 HZ] , whose result is the starting point of this theorem. From the holographic algorithms 
point of view, most of the known holographic algorithms (301 are essentially for planar 2-3 regular graphs. 
The dichotomy theorem here explains the reason why they are special and why many variations of them are 
#P-hard. In the previous two dichotomies for Pl-Holant c and P1-#CSP, the new tractable cases for planar are 
also done by holographic algorithms with matchgates. However, only special basis transformations are used since 
we assume some signatures are freely available. In this planar 2-3 regular graphs setting, no additional signatures 
are assumed to be freely available. Therefore all possible bases can be used in tractable cases. 

Theorem 6.1. Let [yo 5 Z/ij 2/2] an d l x 0: x ii x 2i £3] be two complex symmetric signatures with arity 2 and 3 respec- 
tively. Then Pl-Holant([j/o, 2/i> 2/2] | [xo, x\, X2, £3]) is #P-hard unless [3/012/1,2/2] and [xo, x±, X2, £3] satisfy one of 
the following conditions, in which case it is tractable: 

1. Holant([j/0j 2/i> 2/2] |[^0i x\, X2, 2:3]) is tractable (for which we have an effective dichotomy [9]); or 

2. There exists a basis T such that both [j/o, 2/ij D2\{T~ 1 )® 2 and T® 3 [xo, Xi, X2, X3] are realizable by some match- 
gates (for which we have a complete characterization \1 Vj). 

Proof. If [xo, xi, X2, X3] or [j/q, Hi, 2/2] is degenerate, the problem is tractable, even for the non-planar case, 
and so this falls in condition 1. Now we assume that they are both non-degenerate. As proved in |16j . we 
can choose an invertible Tj such that [xo, Xi, x 2 , X3] (or its reversal, which is similar and we omit that case) 
can be written as T® [1,0, 0,1] or T® [1,1,0, 0]. Therefore by a holographic reduction, we can always reduce 
the problem equivalently to one of the following two problems: (1) Pl-Holant([zo, z\, Zq\\ [1, 0, 0, 1]) and (2) Pl- 
Holant([zo, z\, ^2] | [1, 1,0,0]). So it is sufficient to prove the theorem for these two cases. 

For Pl-HolantQzo, Zi,z%] | [1,0,0, 1]), by the theorem |2.11| ([27]). the only case which is hard for general graphs 
and tractable for planar graphs is zf t — z\. This condition is exactly the same as the condition that there exists 
a basis T such that both [y , yi, 2/2] (7 1-1 ) 182 and T® 3 [1,0, 0, 1] are realizable by some matchgates. This proves 



Theorem 6.1 for case (1). 

Now we consider Pl-HolantQzo, z\, zi\ | [1, 1, 0, 0]). If zq — 0, the problem is trivially tractable even for general 
graphs. This can be seen by a simple counting argument: in a bipartite graph the LHS vertices all have the 
signature [0,^1,^2] and thus at least half the edges must be 1, while the RHS vertices all have the signature 
[1, 1, 0, 0] and thus less than half the edges are 1. This is also the only case where the problem is not #P-hard for 
general graphs when the RHS has [1, 1, 0, 0] by [9]. Now we assume zq ^ 0. Then it is sufficient to prove that either 



\l (Z0Z2 



the problem is #P-hard or there exists a basis transformation T such that [1, 1,0, 0]T® 3 and (T 1 )"°"'[zo,zi 



are realizable by some matchgates. Let T — 
invertible since zq ^ and [zq, z\, Z2] is non-degenerate ( i.e. 



(^i) 2 )/^o 



Note that T is well defined and 
Then we can verify that 



[1,1,0,0]T 



®3 



z (z + 3zi), y/zo(zoz 2 - (zi) 2 ), 0, 0] and (T- v f 2 [z Q , z u z 2 ] - [1, 0, 1]. 



We note that \/ z (z Z2 — (zi) 2 ) ^ 0. If ^/zq{zq + 3z ± ) = 0, then both [^/zq{z + 3zi), y/ 'z Q (z z 2 — (z!) 2 ),0,0] 
and [1,0,1] can be realized by matchgates and the problem for planar graphs is tractable. We denote v = 
Vzo(zo+3zi) _^ q Then the problem is equivalent to (non-bipartite) Pl-Holant([w, 1,0,0]). Now it is sufficient 

yf z a (z z 2 -(z 1 ) 2 ) 

to prove the following claim: 

Claim: Let v ^ be a complex number. Then Pl-Holant([w, 1, 0, 0]) is #P-hard. 

We can realize [v 3 + 3v,v 2 + l,v, 1] by connecting 3 copies of [v, 1,0, 0]'s as illustrated in Figure [9j If we 
can prove that Pl-Holant([w 3 + 3w,w 2 + l,v,l]) is #P-hard, then we are done. In tensor product notation this 
signature is 

[v 3 + 3v,v 2 + l,v,l} 1 





V + 1 


03 


' V- 1 " 


®3\ 




1 


+ 


1 





Then the following reduction chain holds: 
Pl-Holant([w 3 + 3v, v 2 + 1, v, 1]) 



Pl-Holant([l,0,l] 



= T Pl-Holant([u 2 



2v 



3v,v 2 + l,u,l]) 

2]|[1, 0,0,1]), 



2,v 2 ,v 2 - 2v 
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Figure 9: All vertex signatures are [v, 1, 0, 0]. 



where the second step is a holographic reduction using 



1 



1 



This transforms the problem to our 



1 1 

first case where the RHS all have [1, 0, 0, 1]. The only possible exceptional case happens when (v 2 + 2v + 2) 3 = 
(v 2 -2v + 2) 3 . Since (v 2 + 2v + 2) 3 - (v 2 - 2v + 2) 3 = 4v(3v 4 + I6v 2 + 12) and v ^ 0, we will have proved the 
claim as long as 3i> 4 + 16u 2 + 12 ^ 0. There are four roots for the equation 3v 4 + l&v 2 + 12 = 0, and for these 
four exceptional values of v, we prove it separately as follows. 

In addition to the gadget in Figure [9j we can construct a gadget in Figure 10 with a binary signature 



2,v, 1]. Now it is enough to prove that Pl-Holant([w + 2,v, 1] 



3v,v 2 + l,v, 1]) is #P-hard. Under 




Figure 10: All vertex signatures are [v, 1, 0, 0]. 

v ~\~ 1 v — 1 

the same basis ^ ^ , we will get an equivalent problem Pl-Holant([A, Y, Z] | [1,0, 0, 1]), where X = 

(v 2 + 2)(v 2 + 2v + 1) + 2v(v + 1) + 1, Y = (v 2 + 2)(v 2 - 1) + 2v 2 + 1, and Z = (v 2 + 2){v 2 - 2v + 1) + 2v(v - 1) + 1. 
Again this transforms the problem to our first case, and, it is easy to verify that any root of 3i> 4 + I6v 2 + 12 = 
is not a tractable case here. This completes the proof of the claim and also the proof of the theorem. □ 



Appendix: Some Connections to Statistical Physics 

In this section we describe some background and connections from Statistical Physics. Our discussion is necessarily 
a superficial one, both due to our limited knowledge and limitation on space. The purpose is to illustrate that, 
even at such a superficial level, a strong connection exists, and that our complexity results may shed some light 
on the venerable question from physics: Exactly what "systems" can be solved "exactly" and what "systems" are 
"difficult". 

The Ising model was named after Ernst Ising (52] . Wilhelm Lenz invented this model and gave it to his student 
Ising to work on it. The model consists of a discrete set of variables, called spins, that can be assigned one of two 
values (states). These spins are usually placed on a lattice structure or a graph, and each spin interacts with its 
nearest neighbors. 

Denoting the values each spin i can take as = +1 and —1, the energy (the Hamiltonian) of the Ising 
model is H(a) = — X)edge{ij} Ji,j a i a j- The interaction between spins i and j is called ferromagnetic if J^.j > 0, 
antiferromagnetic if Jj j < 0, and noninteracting if J^j = 0. E.g., if all the spins are placed on a one-dimensional 
lattice, then the antiferromagnetic one-dimensional Ising model (with the same value Jjj = J < 0) has the energy 
function H — (7^+1, after normalization. The ferromagnetic two-dimensional Ising model on a square lattice 
(with the same value Jij = J > 0) has energy H = — j( a i,j a i,j+l + The Ising model may be 

modified by magnetic fields which amounts to a unary function at each spin H = — Sedge{i j} ^ i J (Ti ' J i ~ ^■i< T i- 

The model is a statistical model. The central premise of statistical physics is that the probability of each 



22 



configuration a is given by the Boltzmann distribution, e H ( cr )/ kT j e ff ( cr )/ feT i where k is Boltzmann constant 
and T is the (absolute) temperature. This focuses attention on the partition function 



Z 



-H(a)/kT 



Note that the exponential e~ H ^ a ^ kT turns this into a sum-of-product functions exactly as we discussed in #CSP. 

In 1925, Ising solved the one-dimensional Ising model [55]. The 2-dimensional square lattice Ising model with 
zero magnetic field was solved by Onsager in 1944 [31]. Onsager announced the formula for the spontaneous 
magnetization for the two-dimensional model in 1949 but did not give a derivation [32 . C.N.Yang (1952) gave 
the first published proof of this formula |43j . using a limit formula for Fredholm determinants, proved in 1951 by 
Szego in direct response to Onsager's work. There are many extensions to the basic Ising model [301 IT]. 

Another landmark achievement is the exact computation of the number of perfect matchings (dimer problem) 
on any planar graph using Pfaffians. This was independently discovered by Kasteleyn and by Fisher and Temper- 
ley [33, 24, 25 . This problem can also be nicely expressed by a partition function in our Holant framework; where 
this time the Boolean variables are the edges (to include an edge or not), and the local constraint function at 
each vertex is the Exact-One function. Freedman, Lovasz and Schrijver [5T] recently proved that this partition 
function cannot be expressed as a graph homomorphism function, where the vertices are variables as in the Ising 
model. However in the framework of Holant problems we can find a unity for all these problems. 

We note the following. In the paper |11| we gave a complete characterization of matchgate realizable symmetric 
signatures. The following lemma is proved [llj : 

Lemma 6.2. The set of bases under which the signature [xq, X\, x 2 ] is realizable as a signature by some matchgate 
is 



e GL 2 ( 



- 2xipini -) 
or xqPqPx 



x 2 n\ 



0, x pl - 2xip «-o + 
zi(rc pi + n iPo) + x 2 n n 1 



x 2 nl 
= 



Q 



This has the consequence that under the basis 



the signature [x, y, x] is 



7m\= [©■(-. 

realizable by a matchgate, for all values x and y. In terms of the Ising model, when two interacting spins i and 
j take the same assignment value o~i = o~j = ±1, the contribution to the Hamiltonian is —J%j, and when they 
take the opposite assignment o~i — —o~j — ±1, the contribution is Jij. Translating this to the contributions to 
the partition function we get exactly the local constraint evaluation x — e Ji '^ kT when inputs are 00 or 11, and 
y = e~ J * j/ fcT when inputs are 01 and 10. 

Then, the theory of Holographic Algorithms tells us that for planar graphs, this Ising model is exactly solvable 
by a holographic reduction to the FKT algorithm. 

The present paper, especially Theorem |5.1[ tells us why this is exactly where physicists stopped, and attempts 
to generalize this to non-planar systems have not been successful in the past 85 years. 

Sorin Istrail |23j showed that computing the free energy of an arbitrary subgraph of an Ising model on a lattice 
of dimension three or more is NP-hard; see a nice article by Barry Cipra in the SIAM News [17]. A very partial 
list of a great deal of research on this and related models, from a computational complexity perspective, can be 

foundin [a ei m in cna eh cm on on ei ei nansiiii]. 
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